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Abstract 


This paper developes a simplified continuum (continuous with respect to time, stress, etc.) fatigue damage 
model for use in critical design, Life Extending Control and fault prognosis. The work is based on the local strain 
cyclic damage modeling method. New nonlinear explicit equation forms of cyclic damage in terms of stress ampli- 
tude are derived to facilitate the continuum modelling. Stress based continuum models are derived. Extension to 
plastic strain — strain rate models is also presented. Progress toward a non-zero mean stress based is presented. Also 
new nonlinear explicit equation forms in terms of stress amplitude are derived for this case. Application of the 
various models to design, control, and fault prognosis is considered. 


I. Introduction 


A strong motivation for the current activities to develop continuum fatigue damage models comes from the 
ongoing work in the areas of life extending and or damage mitigating controls (refs. 1 to 3, respectively). These 
controls studies seek to create control methodologies to allow the reduction of damage in critical components in 
aerospace systems by the manner in which the control moves the system transiently between setpoints. The transient 
damage for critical components in rocket engines has been shown to be capable of reductions on the order of 2/3, by 
the manner in which the control moves the engine through the transient. This has been accomplished without 
significant loss in dynamic response (ref. 4). The results quoted above are based on open-loop studies which have 
been accomplished through the use of gross nonlinear optimization. The continued development of life-extending 
control requires a damage model which is continuum based, as opposed to current fatigue damage models which are 
cyclic extrema based. Contemporary cyclic methods require a completed stress strain cycle before the associated 
damage can be determined. What is required for controls is the ability to predict for the next increment, or con- 
tinuum, of stress or strain what the associated fatigue damage will be. That is the thrust of this paper. Clearly it is 
desirable to create as simple a continuum model as possible, since this will allow a more broad application of the life 
extending control concepts. Only very limited work has been done in this area. The initial work done by A. Ray et al. 
at Penn State University (ref. 5) has created a useable continuum model for fatigue damage. However, this model is 
hindered by the requirement of identifying the cycle extrema and then calculating the damage between the extrema 
based on extrema information. This complicates considerably the use of such a model in a practical control design. It 
requires a continuous accounting of the cycle extrema and the changing (bookkeeping) of these extrema as the 
physics of the process progresses. An approach which would use only local stress or strain to infer damage would be 
simpler and hence superior. 

This paper seeks to create a zero mean stress continuum fatigue damage model in two forms, the first form is 
stress based and the second is plastic strain and strain rate based. The paper introduces new simplified forms for the 
cyclic damage results, for the zero mean stress case. Progress on a non zero mean stress continuum model is also 
shown and the open issues that remain in this area are discussed. 

The basic objective of this work is to generate a damage model for the fatigue failure of metallic materials 
which is continuum or differential based as opposed to current theory which is cycle based. The fundamental ap- 
proach of the paper is to use results from the cyclic local strain method as a basis for the development of the con- 
tinuum model. The next section will present a short summary on those parts of the local strain approach required for 
the material which follows. The work of Dowling et al. (ref. 6) will be the focal point that will be used in this paper 
as the basis of the local strain method. This paper expands on the work in reference 7 in particular with regard to 
applications. 
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Local Strain Method 


The basis of current damage (fatigue/fracture) approaches, study the experimental results of applying cyclic 
loads of various amplitude and various mean (constant) bias loads, and summarize/generalize these to allow predic- 
tion of arbitrary combinations of loading cycles. Many variations and methods have been evolved. A fairly straight 
forward approach (called the local strain approach) by Dowling et al. (ref. 6) will be the foundation for the analysis 
which follows. 

Atypical stress-strain hysteresis loop is shown in figure 1(a). The effect of cycle amplitude changes (with 
zero mean stress) is illustrated in figure 1(b). The back bone plot (cyclic stress-strain curve) which is the locus of the 
extrema of the stress-strain cycles is shown (ref. 6) to have the mathematical form (for materials of interest); 
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Figure 1. — Cyclic stress-strain, (a) Typical hysteresis loop, (b) Stress-strain 
cycles of various amplitudes-with backbone curve. 
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where e a =Ae/2 and G fl = Ag/ 2 are strain and stress amplitude, respectively, and£. A, and s are material 
constants. The cyclic damage associated with repeated hysteresis stress-strain cycles of a given amplitude is deter- 
mined by experimental observation to have the mathematical form: 

e «=^K)‘ +e /Kr <2) 


where b , c, E , Oy- and are constants for a particular material and N j is the number of cycles to failure. The plot 
of figure 2 shows a typical curve as described by equation (2). 

The total strain amplitude (eq. (1)) is seen to be composed of two parts the elastic strain contribution. 


and the plastic strain contribution 
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Equation (2) now can be used to estimate the damage associated with a cycle of strain of amplitude e fl . Since Nf is 
the number of cycles to failure at amplitude e a then 1/ Nj is the damage of a single cycle (assuming no effect of 
accumulated damage). Therefore n cycles of amplitude G a will create a damage D of 



Figure 2. — Cycles to failure versus strain amplitude. 
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where D = 1 represents failure. Then the Palmgren-Miner equation may be used to determine the damage for cycles 
of different amplitudes, i.e. 



( 6 ) 


where i represents the various amplitudes composing the strain history. Various corrections can be applied for the 
effect of mean stress superimposed on the cycles and alternate methods have been evolved to account for the effect 
of damage accumulation (nonlinearity) on the damage of any particular cycle. 


Zero Mean Stress Continuum Damage Model 

Because of the strong nonlinear terms in equations (1) and (2) above, these equation forms are not directly 
suitable for the continuum model development. Fortunately, there is a redundancy in the material properties 
( b, c t E, c'p z' f , and s) that can be used. This can be determined from equations (3) and (4) by solving for 2 Nj in 

each and eliminating 2 Np Thus 


e =e> 
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\c/b 
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where £ is the amplitude of the plastic strain. 
Pa 

Now using (a a /A) 1/j to eliminate z p 



and rewriting the right hand side of equation (7) gives 
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For this to be true the following must apply 


1 _ c 
s b 


( 9 ) 


and 


Ae' f b/C = Ae'f = o' f (10) 

For the materials study in reference 6, these relationships are found to be virtually exact for RQC-100 and approxi- 
mately correct for MAN-TEN steel. These equations are assumed to be generally approximately correct for many 
metallic materials and will be used in the analysis which follows. 
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In the analysis that follows stress, o\ will be considered as the independent variable. While it is true that in 
typical fatigue testing the strain amplitude is held constant, it also seems clear by analogy to fluid and current flow 
and other physical processes that load or potential (stress) is the cause of motion (strain in this case). Therefore, 
equations (1) and (2) will be combined eliminating the strain amplitude e fl , thus 
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or since 5 = 1 / N f , where 5 is the damage per cycle this can be written as 
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The elastic terms and the plastic terms in equation (12) are now considered separately to determine 5 cyc . For the 
elastic terms, 
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Using the material relations of equations (9) and (10), the damage per cycle for the plastic part of equation (12) 
becomes 
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Figure 3. — Damage per cycle, 8 cyc versus stress amplitude <r a , for Manten and 
RQC-100 steels; for base computation (EQ 12) and explicit computation (EQ 15). 


Now since equations (13b) and (14c) are the same and satisfy both parts of equation (12), it is the solution (restate- 
ment) of equation (12). Thus, the zero mean stress case cyclic damage law can be written. 
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These equations represent the average (or midlife) damage/cycle in terms of the stress amplitude (at midlife) for 
strain controlled cyclic fatigue. To validate this result against data, the versus Nf life results for RQC-100 and 
MANTEN steels (ref. 6) were used as a test case. To do this equation (12) was solved numerically for 
o a versus 8 and compared to the results of equation (15), figure 3. The results appear to be within the accuracy 

of the numerical solutions over five decades of 8 . The explicit form of equation (15) now allows further analysis 

toward a continuum damage model. 

With stress as the independent variable driving damage and knowing the damage per cycle (eq. (15)) then 


J 8'(a)rfa = 8 cyc 

cycle 
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where 8'(a) is the damage rate 8'(a) = dh! da as the cycle is trans versed. Now a critical question is; over which 
part of the cycle does the damage occur. Various assumptions are possible, with reference to figure 4. It is plausible 
that damage is not likely generated during the relaxation (unloading) phases of the cycle, namely A — > B and C — » 
D, although cracks and acoustic emissions may be observed. Damage is most likely during tensile stressing D — > A 
and may also occur during compressive stressing B — > C. 
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Figure 4.— Assumed damage producing domains. 


The initial assumption will be that damage occurs only for cr > 0 and a increasing (domain D -> A, fig. 4), 
again assuming linear damage accummulation then 
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It is readily shown that 
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also works. Thus the rate at which damage accumulates over the cycle is given by equation (18). This result is 
extremely nonlinear, for example for RQC-100 steel, b = 0.075, = 1.68 x 10 5 is 


S'(a) = 1.5872 xlO -4 


n.12.3333 


V1.68X10 3 ) 


(18d) 


This form shows the extreme sensitivity to small changes (and errors) in stress. It is important to note that equation 
(18) is independent of cycle amplitude, i.e., applies to any zero mean stress cycle. Also equation (18) does not 
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Figure 5. — Damage rate versus stress, (a) RQC-100 steel, 
(b) Location on hysteresis cycle. 


depend on knowledge of the cycle extrema (reference values). A plot of equation (18) for RQC-100 steel is pre- 
sented in figure 5(a). Figure 5(b) shows the damage rate location in a cycle. 

If damage occurs equally in tension and compression then 
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applies over domain B — > C and D — » A in figure 4. In view of the mean stress effect on damage, it is more likely 
that greater damage occurs, over the domain D — > A than B — > C, therefore an unequal damage distribution such as 
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Figure 6. — Damage rate distribution for combined tensile and 
compressive damage. 
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may be a more likely scenario than the equal damage case equation (19). In this equation the parameter k , weights 
the tension side damage relative to the compressive side. Figure 6 illustrates the damage distribution for this case. 

Any of the above forms (eqs. (18) to (20)), can be used as a basis for damage estimation for Life Extending 
or Damage Mitigating Control. These forms are particularly useful because they depend only on stress (and material 
constants), which can be estimated from associated structural models. Clearly the case represented by equation (20), 
is the most general form containing the other two by proper selection of k. However the case of equation (18) shows 
the simplicity of the approach and is easily used to derive further results which may then be generalized. 

To convert the above results to a model useable in the time domain, the case with damage occurring only 
with increasing tensile load (eq. (18)) will be considered. 

For this case the damage rate is given by: 
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In the time domain the damage rate D will be given as 
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and the accumulated damage will be 
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Extension to the more general case equation (20) is obvious. 

It is also noted that D(t) is a monotonically increasing function of time. 

Strain Strain-Rate Continuum Damage Model 

The Damage Equation (eq. (18)) in terms of stress is based on an analysis of a hysteresis cycle, under the 
assumptions of zero mean stress and that all the damage occurs during the extensive (tensile) part of the loop (a > 0 
and a increasing). Extension to compressive damage is obvious and leads to symmetric terms with fractional 
multipliers (eq. (20)). 

The results are believed to be correct so long as hysteresis cycles (or similar load) are being analyzed. 
However, consider, the loading case shown in figure 7. In this loading scenario the last segment, A — > B, completes 
a hysteresis loop B y F, F, A, B . But shortly into the return (Point C) the loading is again increased. A profile similar 
to C — > D will be experienced. Note that most of C — > D is elastic. After point D is reached, strong plastic strain will 
be again be experienced. 
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Figure 8. — Reference hysteresis cycle for strain/strain rate model. 


Integration of equation (18) will predict the same damage for C — » D as A — » B, since the stress levels and 
changes are the same. Use of equation (18) for A — » B will yield a correct damage estimate. However, it is believed 
that its use on leg C — » D will be overly conservative since the plastic deformation is relatively much smaller. It is 
further believed that a damage rate equation based on plastic strain and strain rate will more accurately predict the 
damage rates on both A — » B and C -» D, and will be more generally applicable. 

The derivation of the damage model in terms of plastic strain and strain rate will be derived for the tensile 
damage only case and is based therefore on equation (18a). 

Consider a hysteresis loop (fig. 8) of amplitude g q1 z a . For the case described only the lower curve FGC 
need be considered. The equation for this curve, reference 6 is given as 


£ ~ £ r 
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where o and e are the stress and strain respectively on FGC and (o r , e r ) are the coordinates of the previous strain 
reversal (point F here). The curve for the locus of reversal points (for (-,-) quadrant) is 
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Letting £ r = - z a and G r = - o a in equation (24) and replacing z a by equation (25) gives 
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where z p is the plastic strain component and z e is the elastic strain component (a/E). Considering only the plastic 
strain component gives 


e p = 2 
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Now, differentiating with respect to g yields the plastic strain rate as 
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Solving for o a from this equation gives 


c a =2A(ASz' p ) -G (29) 

This equation indicates that given a family of zero mean stress hysteresis curves, of various amplitudes, the particu- 
lar curve corresponding to stress amplitude, o a , can be identified knowing the plastic strain slope z' p and the stress 
level. 


Now o a can be replaced in equation (27) by the expression of equation (29) and solving for a gives 
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This equation relates the stress level to the plastic strain rate and strain over curve FGC. The relationship holds 
regardless of cycle stress amplitude o a . The damage rate 5'(o) in equation (18) can now be expressed in terms of 

z p and e^by substituting equation (30) into equation (18), thus 
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a > 0, a increasing 


This is the basic result of this section, it applies to the tensile stress only case. Extension to the compressive damage 
cases requires consideration of the curve CHF in figure 8. Again reference 6 gives the equation for this curve as 
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With CT r = CT a and e r = e a = o a /E + (o^/a)*^ the plastic strain for CHF becomes 


fn \ 1/s ( rr -rr\ 1/j 


£ P = 


V A J 


q g~ g 

2A 


( 33 ) 


and 


Ep As 




2A 


for CHF 


(34) 


Following a derivation similar to that above and using equation (19) as the basis gives 
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Figure 9. — Damage computation based on plastic strain and plastic 
strain rate. 


13 






as the compressive damage only expression. Extension to the combined tensile compressive damage case of equa- 
tion merely requires the multiplication of equations (31) and (35) by k and (1 - k) respectively. 

The discussion that follows shows how z' p j can be determined in application. It will be assumed 

stress, a, and strain, e, measurements or estimates are available at the critical load (damage) point(s) of the structure. 
It is desired to estimate z p (t) and z' p (t) at such a point. The estimate(s) are based on the graph of figure 9. The 
Elastic Modulus ( E) is assumed to be known and constant. During elastic straining (loading) 
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Then the plastic strain at point A, is given by 


e 


Pa 



(36) 


The plastic strain rate z' p is determined based on a small step from point B, thus: 
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This form is valid for both legs of the hysteresis cycle so long as B is taken as the trailing point. 

Now of course equations (31) and (35) are converted to time dependent forms as was done in equation (22) 
i.e., 
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Non-Zero Mean Stress Continuum Damage Model 

The previous sections have suggested continuous models based on the case of zero mean stress. It is 
important to examine the effect of mean stress in this regard. While many expressions have been generated in the 
cyclic damage format for the effect of mean stress, the following equation (ref. 8) will be used to attempt to generate 
a continuum nonzero mean stress damage model: 
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In this expression G m represents the mean cyclic stress. Assuming the cyclic stress-strain behavior of the material is 
not altered by mean stress then 
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and Ay = 1/8^ may be used to express equation (39) in terms of stress amplitude and damage per cycle giving 
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Equating the plastic parts of equation (41) gives 
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Using the material relationships, equations (9) and (10) and after some algebraic manipulation, solving for 8 ( 
yields amazingly 
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The logic proceeds as in the zero mean stress case, i.e., since 
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satisfies both elastic and plastic terms of equation (41) it is a solution of equation (41) and represents a simplified 
mean stress cyclic damage law. The value of this form (eq. (46)) over that of equations (41) or (39) is that it may be 
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Figure 10. — Effects of mean stress and stress amplitude on damage 
per cycle, 5^, for RQC-1 00 steel. 



Figure 1 1 . — Assummed damage location for tensile side only mean 
stress damage model. 


explicitly solved for any quantity. Figure 10 shows the character of equation (46) for various levels of damage/cycle, 

5 eye 

With the availability of the explicit form (eq. (46)) a mean stress condition continuum damage model 
formulation may be attempted. 

The case where all the cyclic damage occurs only between a m and G m +a a on the tensile leg (fig. 11) is 
considered. The problem is to determine a damage rate 5'(a) over this domain such that: 
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the problem can be reformulated as 


J o ° fl 5'(o+a m ) Ja = 2 
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Then differentiating the right-hand side of equation (48) gives 
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and the desired integrand then is given by 
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This then (eq. (49)), is the desired damage rate for a hysteresis cycle with a mean stress condition. The hope prior to 
starting this part of the analysis was that this expression would be found to be independent of o m so that a “univer- 
sal” stress based continuum form would be obtained, this unfortunately was not the case. 

This equation (eq. (49)) parallels equation (18) and can be seen to reduce to it when o m = 0 . Further, if it 
assumed that damage occurs equally in tension increasing and compression increasing modes, then 
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may be used, paralleling equation (19). And for unequal damage on the tensile and compressive legs for the mean 
stress case 
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and equation (51) parallels equation (20). 
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(51) 
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Figure 12. — Damage rate verses normalized stress level with 
= 0 , . 25 , . 50 , . 75 . 

The effect of mean stress level on damage rate for RQC-100 steel is shown in figure 12 for the tensile only 
damage assumption (eq. (49)). It is very important to note here that for stress levels less than 0.925 the zero mean 
stress damage rate is greater than all non zero mean stress cases (for a m =0). That is, the simple zero mean stress 
damage law of equation (1 8), is conservative in predicting damage rate at these stress levels (with mean stress). This 
suggests that it may be feasible to use the envelope (maximum) damage curve as a conservative damage law for 
practical Life Extending Control applications. 

For a m > 0, the maximum damage curve at a given (constant) stress level, is determined by setting the 
derivative of 5 '(a) in equation (49) to zero, thus 
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After considerable algebra, this yields the condition 
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Substituting this results into equation (49) and after considerable simplification 
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for a' f > > 0 
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(54) 


In words; the maximum damage rate is inversely proportional to the distance of the stress from ay . It can be shown 
from equation (53), with o m l = 0 that equation (54) only applies for 




(55) 


Thus for RQC-100, l+b = 0.925, and equation (54) is seen to determine 8^ for a/a' f > 0.925 while below that 
level the zero mean stress equation (18) dominates. Extension, of the mean stress based damage laws to the time 
domain follows that shown for equations (22) and (23) with obvious changes in db/da. The plot of damage rate 
versus normalized stress level (fig. 13) compares the mean stress results with the maximum damage rate result 







Figure 13. — Damage rate verses normalized stress level with 
= 0 , . 25 , . 50 , . 75 . 


19 



predicted by equation 54. Note, if G m is allowed to be negative damage rates can exceed the zero mean stress 
damage prediction. 

Applications 

Controls Applications 

The original intention of this work was to provide a simple continuum fatigue damage model for use in 
damage mitigating (for life extending) control. Two areas of application are possible; control of accumulated 
damage thru transients and secondly the avoidance of steady state operation in operating domains of high damage 
(when possible). To be useful for controls the primary requirement for a fatigue damage model is that it is continuum 
based (as opposed to cyclic), beyond this it must properly indicate at least approximately the magnitude of damage 
with respect to some primary variable (here stress). That is to say that the damage model must always indicate the 
correct direction to move the system to reduce damage however, the exact magnitude of the system movement (thru 
high model accuracy) may not be so important. 

The damage models derived here (eqs. (18), (19), and (20)) are simple in that they are dependent only on the 
local stress level. A step more complex are the mean stress based models (eqs. (49), (50), and (51)). Here both stress 
and mean stress must be known to estimate damage rate. The ultimate in simplicity is the maximum damage 
envelopes model (eq. (54)). It depends only on stress and will always yield a conservative result and provide a 
proper direction for Life Extending Control. 

Design Applications 

Availability of simple continuum fatigue models allows the multidisciplinary optimization of hardware 
design and system operation for long life. This concept has been forwarded for Reusable Rocket Engines (ref. 9) 
where it is proposed to optimize both aspects simultaneously. For design and operation optimization, again simplic- 
ity of computation is important since many computer iterations will be necessary. Use of a conservative model such 
as equation (54) could be used to closely approach an optimum then for final refinements more accurate (less 
conservative) model could be substituted to complete the optimization. 

Prognostic Applications 

The determinations of critical component health could be done using the strain-strain rate model (eqs. (31) 
and (35)), it is believed (not verified) that this would provide more accurate results than the stress based models 
however, both stress and strain information would be required and the forms are considerably more complex than 
those that are stress based. The conservative forms (eqs. (20), (51), and (54)) may be used directly and could provide 
conservative prognosis of life/impending failure. Further it may be possible to calibrate these models for dedicated 
application to improve accuracy while maintaining their simplicity. It is noted that for repetitive cyclic damage these 
forms will be as accurate as the local strain method on which they were based. Use of the non-zero mean stress 
damage models (eqs. (49), (50), and (51)) pose the additional burden of determination of the mean stress a m for each 
cycle. For this reason use of the maximum damage rate form (eq. (54)) may be preferred even though it is 
conservative. 

Concluding Remarks 

This effort seeks to create a continuum fatigue damage model. This initial effort attempts to mathematically 
convert the classic cycle based damage results as represented by the Local Strain Approach into a continuum model. 
The results achieved to date have been encouraging. Under the assumption of zero mean stress it has been shown 
that the continuum damage rate can be expressed as a highly nonlinear function of the instantaneous stress (eqs. (18) 
to (22)). These continuum rates are compatible with hysteretic damage. In the process of achieving these results 
explicit simplified forms of the cyclic damage laws have been developed (eqs. (15) and (46)). Also time domain 
implementations have been derived (eqs. (22) and (23)). A continuum fatigue damage model in terms of plastic 
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strain and plastic strain rate has been shown to be equivalent to the stress based model for zero mean stress hyster- 
etic damage but is believed to be superior for general application. This has not been validated. 

The case of non-zero mean stress is more difficult both conceptually and analytically. Here also, an explicit 
simplified form of the cyclic damage law (eq. (46)) has been developed. Based on this explicit form, stress level 
based damage rate expressions (eqs. (49) to (51)) have been derived under various assumptions. The assumptions 
center around broad questions of where in the cycle damage occurs. 

Initial comparisons between the positive mean stress case and the zero mean stress case have shown that for 
stress levels up to 92.5 percent of that the zero mean stress continuum model yields conservative results (i.e., 
greater damage). Thus for Life Extending (or Damage Mitigating) Control it may be applied to most practical 
situations. 

Left open is the question of the nonlinear damage effects, namely the effect of accumulated damage on 
damage rate. While not reported yet, some progress has been made in this area with the availability of the simplified 
expressions discussed above. What is needed are appropriate cyclic data sets showing variations in stress amplitude 
versus cycles, for constant strain amplitude testing, to allow calibration of the effect. Further research effort is also 
needed to create plastic strain/strain-rate damage model for the mean stress case. Further, study of the relationship 
between this model and that for the zero mean stress should be done. 

Finally experimental studies are required to screen broadly where the intracycle damage occurs. Such 
studies are now being started. Availability of such data will determine which assumptions of the analysis (i.e., 
tension damage only, etc.) apply. In the longer run further experimental validation of the proposed continuum 
models will also be required. 
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